Observable geometric phase induced by a cyclically evolving dissipative process 
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In [J a new way to generate an observable geometric phase on a quantum system by means of a 
completely incoherent phenomenon has been proposed. The basic idea is to force the ground state 
of the system to evolve ciclically by " adiabatically" manipulating the environment with which it 
interacts. The specific scheme analyzed in Q], consisting of a multilevel atom interacting with a 
broad-band squeezed vacuum bosonic bath whose squeezing parameters are smoothly changed in 
time along a closed loop, is here solved in a more direct way. This new solution emphasizes how the 
geometric phase on the ground state of the system is indeed due to a purely incoherent dynamics. 
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INTRODUCTION 

Geometric phases appear when a quantum state un- 
dergoes a cyclic evolution. The most celebrated example 
of geometric phase is perhaps the Berry phase 0]. Let 
us illustrate the phenomenon with a simple system: a 
spin in a magnetic field whose direction changes slowly 
in time. The Hamiltonian of such system is 

H=-B(T (1) 

where er = (<t x , a y , a z ) are Pauli operators and B(t) = 
Bq (sin i9 cos tp, sin sin ip, cos "&) is a three dimensional 
vector, which we assume to be time dependent. The en- 
ergy eigestates of H$ are the eigenstates of the operator 
er • n where n = (sin $ cos sin $ sin cost?) is a unit 
vector pointing in the instantaneous B direction and can 
be written as 

I Tn) = cos^| t z ) + e^sin^| Q 
Un)=sin||T,)-e^cos|u,) (2) 

where | i z ),\ tz) are the eigenstates of the operator a z . 
Let us assume that n changes slowly in time so that 
the so called adiabatic approximation is valid and that 
at t = the system is prepared in an energy eingen- 
state, say |^(0))) = | Tn (0))- At at time t we have 
IV'(i))) = I Tn (t)) i.e. the state of the system adiabati- 
cally follows the direction of n. Let us now further as- 
sume that B(T) = B(0) i.e. the adiabatic motion of B 
is cyclical. In this case the energy eigenstrates acquire, 
on top of the dynamical phase, a phase factor Xg that 
depends only on the geometry of the path followed by 
B(t). In mathematical terms 

| Tn (T)) = e<*>e iB ° T \ U (0)) (3) 



where for simplicity h = 1. The geometric phase is given, 
according to Berry [jj, by 

Xg = if(U\d\U) (4) 

A straightforward calculation shows that \g is equal to 
the solid angle spanned by B(t) in his cyclic evolution 
and it is therefore independent on T. Note, by the way, 
that the eigenstates | I) and | f) acquire a geometric 
phase of opposite sign. 

The role of the adiabatic approximation in the example 
above is to provide a way to perform a parallel transport 
of a state vector along a suitable path in Hilbcrt space 
by slowly changing a time dependent system hamilto- 
nian. One may wonder if the same goal can be achieved 
by changing slowly in time the irreversible dynamics of 
a system coupled with a reservoir. To be more specific 
assume that a system interacting with a rigged reservoir 
relaxes to an equilibrium pure state \ip s ). Such state will 
in general depend on the properties of the reservoir. If 
the properties of the reservoir are changed slowly in time 
the equilibrium state \tp s ) follows adiabatically and, if 
such change is cyclical, it acquires a geometric phase Xg ■ 
Furthermore, being induced by the irreversible dynamics 
itself Xg should be intrinsically immune from noise. We 
are not referring here to the standard resistance of geo- 
metric phase to environmental noise which has attracted 
much attention for its potential benefits in the implemen- 
tation of quantum gates §§]. Here something conceptu- 
ally different happens: it is the environment itself which 
generates the the geometric phase. 

A FOUR LEVEL SYSTEM INTERACTING WITH 
A BROADBAND SQUEEZED VACUUM 

The system which we will explicitly analyze here to 
illustrate the above idea is the same studied in 111: a 
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multilevel atom interacting with a broadband squeezed 
vacuum. The interesting feature of such system is that its 
irreversible dynamics relaxes to a nontrivial pure ground 
state, whose configuration depends on the squeezing pa- 
rameters of the reservoir. By changing the squeezing pa- 
rameters it is possible therefore possible to manipulate 
indirectly the ground state of the system. 

For the sake of completeness let us first review the es- 
sential features of the dissipative dynamics of an atomic 
system interacting with a broad band squeezed vacuum 
reservoir. Consider first a three level atom whose inter- 
action with an electromagnetic field in the rotating wave 
approximation is described by the following Hamiltonian: 

H = H S + J a\oj)a{u)u)duj + J [g{u)S ] a{uj) + H.c] duj, 

(5) 

where H$ = Q Ylk=-i k\k)(k\ is the internal Hamiltonian 
of the atom, S = \ — 1)(0| + |0)(1| is the atomic operator 
describing the absorption of an excitation by the atom 
and a(ui) is the annihilation operator of the mode with 
the energy oj. The field, which we treat as a reservoir, is 
assumed to be in a broad band squeezed vacuum state 
defined as 

l va c(f?)) sq = K(r))\v&c), (6) 

where K(r]) is a unitary multimode squeezing trans- 
formation jj, 0, which correlates symmetrical pairs of 
modes around the carrier frequency Q 

K(i]) = exp i [J [r)tf(n - u)a)(n H.c] dcjj , 

(7) 

In the above expression -q = e %v r is the squeezing param- 
eter, whose polar coordinates tp S {0 . . . 2ir} and r > 
are care called phase and amplitude of the squeezing, 
respectively. For a broad band squeezed vacuum, char- 
acterized by a value of squeezing parameters basically 
constant over a broad band, one has the following expec- 
tation value for the field bosonic operators: 

(a 1 {oj)a{uj')) = sinh 2 rS(co - J) (8) 

{a{w)a) V)) = (sinh 2 r + l)5(u - J) (9) 

(a(uj)a(uj')) = (e lv sinh r cosh r)6(u> + J - 2fl) (10) 

(a t (w)a t (w / )) = sinh r cosh r)S(w + - 2Q) (11) 

The use of the Born Markov approximation, justified 
by the broadband nature of the field, leads to a master 
equation for the atomic degrees of freedom which can be 
written in the following compact form 

^ = -^{R^Rp + pR^R-2RpB)} (12) 



whith 

7 = 27r|.g(f!)| 2 (13) 

and where we have defined the new "dressed" atomic 
operator 

Rirf) = S cosh r + e lip S f sinh r. (14) 
From (fH|) it follows that the state 

|^H»?)>=c|-l>-e*»|l>, (15) 
with c(r) = / COS u r and s(r) = * in \ r , satisfies the 

v ' V cosh 2r v ' vcosh 2r 

condition R{r\)\^ DF (j])) = 0. In other words, this state, 
being unaffected by the environment, is decoherence free 
IHUIIl El Moreover \ip(r])} represents the new ground 
state, as all the other states of the atomic system relax 
to it. 

So far we have assumed the squeezing parameter r\ is 
time independent, however one may think of a scenario 
in which r\ is slowly changed in time. In such case the re- 
duced atomic dynamics is described by a time dependent 
master equation: 

J = -\{R\t)R{t)p+ P R\t)R{t)-2R{t)pR\i)\ (16) 

It is reasonable to expect that if such change is made 
slowly enough \ipijp(i](t))) is adiabatically changed ac- 
cordingly and and, if 77(0) = rj(T), acquires a purely ge- 
ometric phase. To show that this is indeed the case we 
consider here, with no loss of generality, the simple case 
in which the squeezing amplitude r is kept constant while 
the squeezing phase ip is slowly changed from to 2n. We 
show that, provided tp is small enough, the state \iPdf) 
is adiabatically decoupled from its orthogonal subspace 
and acquires, after a cyclic evolution, a geometric phase 
depending, in this case, only on the amount of squeezing 
r. Note however that, since the steering process is es- 
sentially incoherent, any phase information acquired by 
a superposition of \iPdf) and a state belonging to the 
orthogonal subspace is inevitably lost, as the latter is 
subject to decoherence. The only way to retain such in- 
formation is to consider an auxiliary level \a), unaffected 
by the noise, playing the role of a reference state in an 
interferometric measurement (see Fig^i) 

For simplicity assume that |o) is unaffected by the en- 
vironment during the whole evolution, and, hence, is time 
independent. The whole information about the geomet- 
rical phase and the coherence retained by the system 
during its evolution is then recorded in the phase and 
amplitude of the density matrix term (a\p\ipDF)- In Q 
the system dynamics was analyzed in a suitable "rotat- 
ing" frame. Although mathematically more elegant such 
approach may, at first glance, mask the fact that the sys- 
tem dynamics is fully incoherent. Here therefore we will 
explicitly solve (|l(jfl with no change of frame of reference. 
We must therefore calculate 
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FIG. 1: Schematic representation of the energy levels of the 
two systems considered. 

(a) A four level systems with an equal energy gap Q between 
states | — 1) and |0) and between |0) and The transi- 
tions between these level are coupled to the modes a(cu) of 
the reservoir. The reference state \a) is decoupled from the 
reservoir. 

(b) A five- level system where transitions 1 <-» and 0^—1 
are coupled to modes 5i(a>) while 1' «-» and <-> —1' are 
coupled to modes &2(uj) of the reservoir. 



^-{a\p\ipDF) = (a\p\ip DF ) + (a\p^\ip D F) 
= -i(pse iv (a\p\l) 
= -i(ps [c{a\p\ip±) - s(a\p\ipDF}] 



(17) 



where we have made use of the fact that R\i(>df) = and 
R\a) = R'\a) — and made the substitution 



with 



e %v [c\ip±) - s\iI>df)] 



\iP±)=s\-l)+e Up c\l). 



(18) 



(19) 



The linear system can be closed with the following equa- 
tion of motion 



^(a\p\ips_) = (a\p\ip±) + (a|p^|V>x) 

= -l(o\pWx) + i<pce i *(a\p\l) 



(20) 



= -gfaMV'-L) + i( P c [c(a\p\ip±] ~ s(a\p\ip DF )] 

where 7 = 7 cosh 2r and we have made use of the fact 
that R iR\ipx) = cosh2r|V^±) and R\a) = Rt\a) = 0. Eq. 
(|17I20|I therefore form the following linear system 



dv 
~dt 



where 



-iG-y 



(a\p\ip DF ) 
(a\p\tp-L) 



G = 



—tps 

(fSC 



ifSC 

-ipe 2 — 17/2 



(21) 



Although such system can be easily solved diagonalis- 
ing G, it is however clear that, in the limit 7 3> ip the 
off-diagonal terms are exponentially suppressed by the 
incoherent term. Hence, the leading contribution to the 
time evolution of the initial \4>df) state is a phase factor 
e lx \ with 



Xg 



2ns 2 



2tt- 



sinh r 
cosh 2r 



(22) 



This demonstrate that, in the lowest order "adiabatic" 
approximation, the system preserves its coherence. In 
fact the amplitude damping of p^ a occurs only when we 
take into account the first order contribution in <p/"f. It 
follows that for small <p the system admits an adiabatic 
limit, in which the subspace H£>p{i) spanned by \^df{^)) 
and I a) is adiabatically decoupled from its orthogonal 
subspace H±(t). For this reason, 7i DJ ?(t) is decoupled 
from the effects of the decoherence, which only affect 
states lying in its orthogonal counterpart. Note that this 
adiabatic limit is different from the " standard" adiabatic 
approximation, where the timescale is set by the energy 
splitting ft between energy levels. Indeed in order to ob- 
serve a geometric phase generated by a slowly varying 
Hamiltonian in the presence of damping one must satisfy 
the condition fl ^> ip ^> 7. In our scheme, as the phase 
is generated by a slowly varying limblandian, one must 
instead satisfy <p -C 7. Within this approximation, then, 
a state prepared in the space H.df(0) is adiabatically 
transported rigidly inside the evolving subspace TLDF(t). 
As a result of this adiabatic steering, when the system 
is brought back to its initial configuration, the coherence 
p^ a acquires a phase \ g that can be interpreted as the ge- 
ometric phase accumulated by the state \tpDF(t))- Note 
indeed that according to the canonical formula for the 
Berry phase, the geometric phase of \tpDFit)) is given by 

Xg = i <j> (lpDF\d\lp D F) = 

f 27T d 
= i I (il>DF\-7-\ipDF)dtp = 2ns 2 . 
Jo dip 

which coincides with 1)22(1 . As anticipated the value of 
Xg depends only on the squeezing, and vanishes as the 
squeezing tends to zero. Moreover, notice that the phase 
Xg is purely geometrical, i.e. there is not any dynamical 
term generated by a Hamiltonian term, as, in absence of 
a streering process, the states inside Rdf have a triv- 
ial dynamics. This feature makes the measurement of 
this phase a relatively easy task. Indeed the absence of 
a dynamics makes unnecessary the use of spin-echo tech- 
niques in an interferometric experiment. 
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A FIVE LEVEL SYSTEM 



can be cast in the following compact form 



To conclude let us re-analyse the detection scheme 
proposed in Q. Let's consider the five- level system as 
shown in picture JJJ. This system consists of two repli- 
cas of the three-level system discussed so far, with the 
level |0) in common. The important ingredient is that 
transitions |0) «-> |1) and | — 1) <-> |0) are coupled to 
a set of field modes different from those inducing tran- 
sitions |0) <-> |1') and | — 1') «-> |0). A simple way to 
achieve this, is to choose, for example, polarisation selec- 
tive transitions, say, left-circularly polarised modes for 
the former transitions and right-circularly polarised for 
the latter ones. The Hamiltonian of such system is 



H = Hs + ^3 / al(u>)ai(oj)uidu> 
+ ^ f [g^Sjdiiuj) + H, 

A — 1 1 " 



(23) 



where H s = Q EiL-i k (\ k ) ( k \ + W) ( fc 'l), and Si = 
| - 1)<0| + |0)(1| and S 2 = | - l')(0| + |0)<1'|, and 
is the annihilation operator of the mode with the energy 
u> and polarisation i — {1,2}. Assume that both set of 
modes Si (a;) and di (u>) are in a broad band squeezed 
vacuum, although with different squeezing parameters 
r/i = r\e lipi and r\i = r 2 e 



«V2 . 



vac 



(Vi,m)) sq = Ki(vi)K 2 (V2)\ 



vac) 



(24) 



where Ki(r)i) are the analogous of the operator J7J) acting 
on modes dj. Under the same assumptions described in 
the previous section we obtain a master equation which 



J 



dp 



(25) 



where 



and 



7i = 2 7 r|< 7i (n)| 2 



Ri{r]i) = Si coshr 4 + e lVi S] smhn 

Such system admits a two-dimensional decoherent-frec 
subspace, spanned by states \ip\) and 1-02) defined in the 
analogous way as \tPdf) of equation (JT3J, where c and s 
are replaced by the corresponding c; = cosh / \J cosh 2r i 
and Si = sinh r, / Vcosh2ri, i — 1,2 . 

Again we assume time dependent squeezing parame- 
ters r\\. In particular let us assume constant r\ ^ r 2 and 

01 = 02- 

Suppose that the system is initially prepared in a co- 
herent superposition of state |V ; i)( 7 ?i) an d 102 (^2)), for 
example: 



|0(O)) 



1 



^(\Mv°i)) + \Mv°2))) 



(26) 



On the light of the previous consideration, it is expected 
that, when the parameters (pi,f 2 close a loop, at time 
t = T = 2ir/(p, the coherence has gained a phase which 
is the difference between the expected geometric phases 
X g i — Xg2 acquired by the states |0i), 1020 respectively. 
To prove this we should calculate the time dependence of 
(0i |/?|02) . To this end, we follow the analogous procedure 
of the previous section, which yields to the equation: 



d(ipi\p\ipi 
dt 



= (lpl\p\lp2) 

= -i<P [(4 



H^l|/#2> + <0l|p|02> 
S|)(^1|P|^2) - Sid(lp^\p\ll}2 



S2C 2 (0l|/o|^)] 



r 



where 00 are defined in analogy to Eq. I|19(l . This equa- 
tion has been obtained by inserting the master equa- 
tion l|24l) and using the following properties: i?i,2|0i,2) = 
0, and RiyLj\ipj-) = 0. To close the equation of motion, 
we need to calculate the derivative of the other matrix 
density elements. Eventually, making use of the fact that 

R UM) = R UMi) = and R l R ^) = cosh2r 4 |00) 
leads to a system of equation, which can be expressed in 
the following compact form 



dv 
~dt 



-iK ■ v 



where 



/ (0l|H02> \ 

(0rH0 2 ) 



and where 
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K 



if x 



/ 4-4 
S2C2 

-SlCl 





S2C2 

-r 2 - 
c 2 



-SiCx 



; 02 

"20 







; 7 2 
20 





-S1C1 

S2C2 



(27) 



S2C2 



: 7i +72 
' 20 



with 7, = cosh 2^7. This apparently complicated ex- 
pression, can be cast in a tensor product structure, much 
easier to interpret. Infact K — 1\ (g> G2 + G\ (8 I2 , where 
Gi are 2x2 matrices analogous to <|21|) : 

Q —( ^4 -<p8\Ci \ 

1 y -^SlCl <y3C^ - 171/2 y ' 

and 

Q = ( ~^4 ^ S 2C 2 \ 

2 V ^ S2C2 ~^ c 2 ~ ^2/2 y ' 

We can thus apply step by step the reasoning of the 
previous section, and infer the evolution in the adiabatic 
limit 7^ 3> if. In this limit, the coherence (tpi \p\1p2) is 
not reduces in module, and acquires a phase \g which is 
the difference of the geometric phase acquired by 4>\ and 

1p2- 



CONCLUSIONS 

In the above sections we have described a novel scheme 
to generate geometric phases by cyclically modifying the 
irreversible dynamics of a quantum system. This can be 
seen as a parallel transport of a decoherence free subspace 
generated by a cyclic change of a rigged reservoir. The 
reservoir we consider is a broadband squeezed vacuum. 
The scheme can be discussed in more general terms, as 
shown in . One can think to extend the above scheme 
to other specific scenario. 

Again we stress that the interesting feature of the 
phases so generated is that they are intrinsically immune 
to noise. Furthermore in the scheme analyzed above a 
specific detection procedure which does not need any spin 
echo has been proposed. 



Xg = Xx- X2 



with Xi = 27T 



sinh 2 r 
cosh 2r 



As in the previous scheme, the visibility is only re- 
duced by a factor which is linear in the "adiabatic pa- 
rameters" tp/^/i, which guarantees the existence of the 
adiabatic limit. The advantage of this modified scheme 
is that the value of the geometric phases can be readily 
measured from the polarisation of the light emitted when 
the system relaxes. Infact, if the value of the squeezing 
parameters are suddenly changed from to zero, the 
states \ipi) are no longer decoherent free, and and decay 
to the the ground states | — 1) and | — 1'). Such dissi- 
pative process is accompanied by two photon emissions 
into the reservoir. Due to the structure of the interac- 
tion (|23|l with the reservoir, the photons emitted due to 
the transitions |1) — > |0) and |1') — > |0), are polarised 
according to the geometric phase accumulated between 
\tpl) and 1^2)- F° r example, if and a 2 {uj) are right 

and left circularly polarised modes, respectively, the first 
emitted photon will be linearly polarized along a polar- 
ization plane depending on Xgi — Xg2- 

|V>i) +e l(Xsl - X92) |^2> -> \R) +e l{x <> 1 - x ^ ) \L). (28) 

The detection of of the emitted photon provides therefore 
a simple detection scheme of the geometric phase. 
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